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We predict surface electromagnetic waves propagating across the layers of intrinsic 
Josephson junctions. We find the spectrum of the surface waves and study the 
distribution of the electromagnetic field inside and outside the superconductor. The 
profile of the amplitude oscillations of the electric field component of such waves is 
peculiar: initially, it increases toward the center of the superconductor and, after 
reaching a crossover point, decreases exponentially. 

PACS numbers: 74.78.Fk, 74.50. +r, 74.72.Hs 

I. INTRODUCTION 

A very recent burst of interest to layered high-T c is due to the discovery of a new 
generation of superconductors based on FeAs layers, LaOi-xFxFeAs 1 - and other such sys- 
tems which possess a similar structure. The conventional representative of layered high-T c , 
Bi2Sr2CaCu20s+x (Bi2212) and Ti2Ba2CaCu208+x (T12212) have a structure of supercon- 
ducting CuC>2 layers with Josephson coupling between them. The layered structures of 
high-T c favor the propagation of so-called Josephson plasma waves^, propagating with 
frequencies above the Josephson plasma frequency uj p . The gap structure of the Joseph- 
son plasma excitation spectra has been experimentally observed from measurements of the 
Josephson plasma resonance^. Josephson plasma waves can exhibit remarkable features, 
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including the slowing down of light, self-focusing effects^ and are linked to applications in 
the THz frequency ranged 

It was recently predicted that the layered structure of high-T c superconductors allows the 
propagation of surface waves^. Such waves propagate below the Josephson plasma frequency 
uj p and propagate in the vicinity of the superconducting surface along the layers. In this 
paper we show that there exist surface electromagnetic TM-waves propagating across the 
superconducting layers. 

The electric, E = {E x , 0, E z }, and magnetic, H = {0,H, 0}, components of the electro- 
magnetic waves are proportional to exp[i (qx — ut)} and decay both in the vacuum and inside 
the layered superconductor. Such surface waves across the layers are strongly influenced by 
an external magnetic field h applied along the superconducting layers. Here we describe 
the propagation across the layers of such surface waves and estimate the influence of an 
external magnetic field on their spectrum. 

Z 
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FIG. 1: Interface between vacuum (z > 0) and a layered superconductor {z < 0) in an external 
magnetic field ho. A layered high-T c superconductor has a structure of superconducting layers 
coupled via intrinsic Josephson junctions. 

II. MODEL AND RESULTS 

Consider an interface between the vacuum (z > in Fig. [T|) and a layered superconductor 
[z < 0). Let the c-axis of the superconductor be along the x-axis so that the vacuum- 
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superconductor interface lies in the :q/-plane and an external magnetic field h is applied 
along the y-axis, parallel to the superconducting layers, see Fig. [TJ The electromagnetic field 
inside the layered superconductor [z < 0) is determined by the distribution of the gauge 
invariant phase difference (p(x,z,t) of the order parameter between neighboring layers. It 
is described by a set of coupled sine-Gordon equations^, that in the continuum limit (see, 
e.g.,-) can be written as, 

(*-*£)(£+'H- A «£- a (1) 

Here we neglect the relaxation terms caused by the quasiparticle conductivity; \ a b and 
A c = c/ojpyfe are the magnetic penetration depths across and along layers, respectively, 
uj p = (SneDjc/he) 1 / 2 is the Josephson plasma frequency. The latter is determined by the 
critical Josephson current j c , the interlayer dielectric constant e, and the spatial period of the 
layered structure D. The gradient of the superconducting phase is related to the magnetic 
field h(z), directed along y, as (e.g.,-) 

32 



dip _ 2ttD_ ( 2 d- 



h{z) (2) 



Using d2J) as a boundary condition at z = 0, we solve Eq. flTJ to obtain the dependence of 
the superconducting phase ^p{z) on the distance z from the interface under a homogeneous 
stationary magnetic field ho, 



(po(z) = — 4arctan 



exp 



~o 



A 

where a positive constant Zq > is defined by the boundary condition 

d<p (z). 2nD 

a 2=0 — ~a — n ° 

oz <P 

so that 



z < 0, (3) 



^0 / 7lD\ c 

Here we study the case of relatively small fields, when ho is less than the critical value h c 
and Josephson vortices do not penetrate the superconductor. 
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A. Surface waves at ho < h c 

We take into account the t and x dependence of superconducting phase tp(x, z, t) as small 
variations around the stationary configuration ipo(z) given by Eq. Assuming 

<p(x, z, t) = <p (z) + <f w (x, z, t), 

as a sum of the static and wave terms, we linearize the Eq. (JT]) to obtain 

Substituting ip w (x,z,t) = exp(igx — iut) we derive an ordinary differential equation 
for f (z), 

+ (i _ ft 2 - ^ £ = (4) 

(l + q *\l b )dz* + { cosh 2 [(,-z )/A c ]r ' U 

where we have introduced Q = uj/uj p . Here we are interested in a solution decaying inside 
the layered superconductor: £(2) — > at z — > —00. The equation (BJ has the form of a ID 
Schrodinger equation for a particle with energy 

E(£l) = tt 2 - 1 

in a potential 

U(z) = 



cosh 2 [(z - Zq)/X c ] 

The bound states corresponding to the waves decaying at z — > — 00, can exist for negative 
energies E(Q) < 0, i.e. for Q < 1. One can write an exact solution of Eq. (j3J) in terms of 
the Hypergeometric function, 

£(*) = (1 - ((zff 2 F^-s,e + s + l,e+l, ±±|^) 

where 

£(z) = tanh 

and 



A, 
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s = o (- 1 + \/ 1 + 8( 1 + ^L 



l-n 2 )(l + g 2 A 2 



Q 6 ' 
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We have studied the behavior of the spectrum of surface Josephson plasma waves by means 
of the WKB approximation valid for 



If the inequalities 



Q = q^ah > 1- 



< (l-Q 2 ) <^ 
hi 



are satisfied, there exists a classical turning point z — Zt- According to Eq. 
defined by the equation E{Vt) = U(z t ) that leads to 

1-Q 2 



this point is 



cosh [(z t - z )/X c 



The "wavefunction" £(z) oscillates in the region Zt < z < and exponentially decays at 
— oo < z < z t . After the procedure of matching the "wavef unctions" at the turning point 
by the connecting formulas known from quantum mechanics, we obtain the quasiclassical 
expression for £(z). For the classically-allowed region Zt < z < 0, we have 



A 



cos 



[E(Q)-U(z)] 1/4 
and the underbarrier "wavefunction" for — oo < z < zt, 

A/2 



^\ +g2 [ Z dzWE(n)-u( z >)- 

A c J Zt 



(5) 



p{z) - E(a)] 



1/4 



exp 



l + Q 2 
A. 



dz'^U(z) - E(Sl) 



The waves ^(z)exp(iqx — iut), corresponding to the solution (jSJ), are Josephson plasma 
waves running along the interface of the layered superconductor and across its layers. From 
Eq. (j2J) we obtain the relation of £(z) to amplitudes of the electromagnetic field components 
in the layered superconductor. For the magnetic field, distribution inside the sample, we 
obtain: 



h c X c 



From the ac Josephson relation, Maxwell equations, and substituting A c = c/uj p ^/e and 
h c = Qo/7rD\ c , we obtain the amplitudes of the electric field components: 

$ 



EJz) 



2ttcD 



h c tt 



2v^ 
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and 



Using the Maxwell equations in vacuum, we obtain the dependence of the electromagnetic 
field components outside the superconductor. This gives an exponential decay for positive 

# vac , E™ c , E™ c oc exp{iqx -iut-k v z), z>0 




ky — "T \ / Q 2 



with the decay constant k v = \/ q 2 — uj 2 /c 2 > for q > lu/c. In the WKB regime, when 
Q 1, we obtain 

' „ Q 

X 2 £ ^ab 

as X ab /X c e <C 1. Because of the large X c /X ab 1 and Q > 1, the surface wave decays very 
quickly in vacuum, on the scale ~ X ab /Q, which is much smaller than A c . 

The ratio of amplitudes for the tangential electric and magnetic fields at the interface 
z = +0, above the surface of superconductor, is 

|^- = !% = -vV-^7c 2 . (7) 

In order to derive the dispersion relation for surface Josephson plasma waves, we calculate 
the ratio E x (0)/H(Q) in the superconductor using Eqs. (jSJ) and (jHJ) and then equate this 
ratio to the vacuum impedance Eq. ((7j). This gives 

x ab n 2 Ji + W 



tan 



l + Q< 



f° TV 

/ dz'y/E{Q) - U{z>) + - 

J Zt ^ 



with Q = qX ab . Because X ab /X c e <C 1, this relation can be simplified disregarding the vacuum 
contribution. Thus, 

A//1 A + Q2 dz'^E(Q) - U(z') = 7T (n + ^ , n = l,2,3.... (8) 

A set of dispersion curves for n = 1, 10 is shown in Fig. [2] for two different values of the 
external magnetic field, h /h c = 0.5 (Fig. [2b.) and ho/h c = 0.9 (Fig. [2b). 

The dispersion relation Eq. (jSJ) corresponds to surface waves of an unusual nature. The 
electromagnetic field does not decrease monotonically into the superconductor. Instead, 
the number of oscillations of £(z) with increasing amplitude occur before the exponential 
decrease. An example of the oscillating field E x (z) distribution in surface Josephson plasma 
waves with for parameters h /h c = 0.5, n — 10 and Q = 50 is shown in Fig. [3j 
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n(Q) 




FIG. 2: (a) Dispersion curves for surface Josephson plasma waves at an external magnetic field 
ho/h c = 0.5. The branches correspond to n = 1,2,..., 10 from bottom to top. The spectrum is 
limited from below by the value Q, > (1 — 2/iq//^) 1 / 2 ~ 0.71 (b) Dispersion curves for surface 
Josephson plasma waves at an external magnetic field ho/h c = 0.9. The branches correspond to 
n = 1, 2, 10, from left to right. 

B. Surface waves at small h 

The equation ([8]) does not describe all branches of the spectrum of surface Josephson 
plasma waves. For relatively small magnetic fields and small frequency Q, when the inequal- 
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FIG. 3: Spatial distribution of the field E x (z) of surface Josephson plasma waves calculated for 
ho/h c = 0.5, n = 10, Q = 50. The dashed curve corresponds to the classical turning point 
z = zt, where the WKB approximation fails, z < corresponds to a plasma wave in the layered 
superconductor, while z > corresponds to vacuum. As seen from the figure, the amplitude of 
the surface wave decays very quickly in vacuum on a scale much smaller than A c . Notice that the 
number of nodes of the distribution E x {z) corresponds to the integer parameter n in Eq. J5J). 




FIG. 4: Spatial distribution of the field E x (z) of surface Josephson plasma waves calculated for 
ho = and Q = 0.001. z < corresponds to a plasma wave in the layered superconductor, while 
z > corresponds to the vacuum. 
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0.0005 



Q 
0.001 



FIG. 5: Dispersion curve for surface Josephson plasma waves at different values of an exter- 
nal magnetic field given by Eq. (J9j). The branches correspond to magnetic field values ho/h c = 
0,0.1,0.2, ...,0.7 numbered from top to bottom. 
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E(n) = n 2 -i< 



2h% 
hi 



is fulfilled, the whole superconducting area z < is classically disallowed. In this case, the 
quasiclassical wave function exponentially decreases starting from the boundary, 



£(*) 



A 



\u{z) - Em 



1/1 



exp 



A, 



dz'y/U(z') - E(Q) 



A typical distribution of the wave amplitude for Q = 0.001 and h — is shown in Fig. HJ 
The corresponding dispersion relation can be written as 



l Q 2 - 



^ \ c e^/-E(Sl) - 2h 2 /h 2 c 

This dispersion curve is shown in Fig. [5] for several values of the magnetic field, ho/h, 



(9) 



0, 0.1, 0.7. We have used the parameters: A c /A a h = 500, and e = 20. Note that Eq. (jHJ) is 
valid not only in the WKB approximation and applicable even in the absence of an external 
magnetic field. 
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III. CONCLUSION 

We have predicted Josephson plasma waves propagating along the surface of anisotropic 
high-T c and across its superconducting layers. There exist different modes of such waves 
ordered by the number of nodes n of the amplitude of the electric field component inside 
the superconductor. The profile of the surface waves is unusual: first, the amplitude of the 
oscillations increases inside the superconductor and, after reaching the last node, decreases 
exponentially. 
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